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ABSTRACT 

Weak emission line quasars are a rare and puzzling group of objects. In this 
paper we present one more object of this class found in the Sloan Digital Sky Survey 
(SDSS). The quasar SDSS J094533.99+100950.1, lying at z = 1.66, has practically 
no C IV emission line, a red continuum very similar to the second steepest of the 
quasar composite spectra of Richards et al., is not strongly affected by absorption 
and the Mgll line, although relatively weak, is strong enough to measure the black 
hole mass. The Eddington ratio in this object is about 0.45, and the line properties 
are not consistent with the trends expected at high accretion rates. We propose that 
the most probable explanation of the line properties in this object, and perhaps in all 
weak emission line quasars, is that the quasar activity has just started. A disk wind 
is freshly launched so the low ionization lines which form close to the disk surface are 
already observed but the wind has not yet reached the regions where high ionization 
lines or narrow line components are formed. The relatively high occurrence of such a 
phenomenon may additionally indicate that the quasar active phase consists of several 
sub-phases, each starting with a fresh build-up of the Broad Line Region. 

Key words: galaxies:active - galaxies:quasars:emission lines - galax- 
ies:quasars:absorption lines - galaxies:quasars:individual:SDSS J094533. 99+100950.1 



1 INTRODUCTION 

Strong, broad emission lines are the characteristic signature 
of Active Galactic Nuclei (AGN). However, the Sloan Digi- 
tal Sky Survey (SDSS) finds, among thousands of new AGN, 
rare objects with very weak or almost absent emission lines 
(Fan et al. 1999, Shemmer et al. 2009, Diamond-Stanic et 
al. 2009, Plotkin et al. 2010). Most of these objects, forming 
a new class of weak line quasars (hereafter referred to as 
WLQs), have been found primary at high redshifts (z>2). 
A few sources lying relatively close by have been studied in 
greater detail: PHL 1811 {z = 0.19; Leighly et al. 2007), PG 
1407-1-265 (z = 0.94; McDoweh et al. 1995), and HE 0141- 
3932 (z=1.80; Reimers et al. 2005). The study of WLQs 
optical/UV continuum properties. X-ray and radio emis- 
sion shows no obvious difference from typical emission line 
quasars and no apparent similarity to the weak-line BL Lac 
objects (Shemmer et al. 2009). 

No generally accepted explanation for the weak- 
ness/absence of emission lines has been found so far. Possi- 
bilities so far considered include: 1) dust obscuration or ex- 
treme broad absorption line (BAL) effect, which is unlikely 



since there are no signs of the deep and broad absorption 
troughs in the spectra (Anderson et al. 2001, CoUinge et 
al. 2005) and the X-ray absorbing column is low compared 
with BAL QSOs (Nn < 5x 10^^ cm'^, Shemmer et al. 2009). 
2) Gravitational lensing does not remove the emission lines 
(Shemmer et al. 2006), and strong CHI] and CIV lines are 
found in gravitationally amplified quasars (e.g. Dobrzycki 
et al. 1999, Wayth et al. 2005). 3) Relativistic beaming is a 
good explanation for the absence of strong lines in BL Lac 
objects but WLQs, in contrast to BL Lacs, have bluer opti- 
cal/UV continua, are radio quiet and show no variability or 
strong polarization (Shemmer et al. 2009, Diamond-Stanic 
et al. 2009). 4) The difference in the ionizing continuum 
also faces difficulties since at least the observed part of the 
continuum in those sources looks typical for a quasar. How- 
ever, those objects may still emit more vigorously in the 
unobserved UV, thus having a higher accretion rate than 
the average quasars, which may affect the lines, as e.g. seen 
through the Baldwin effect. No black hole mass determina- 
tion or Eddington ratio was possible for the sample of ob- 
jects analyzed by Shemmer et al. (2009) so as to confirm or 
reject this hypothesis. In this paper we present a new WLQ 
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object found in the SDSS survey which opens a possibihty 
to address this issue. 



2 SELECTION PROCEDURES AND 
OBSERVED PROPERTIES OF THE 
QUASAR SDSS J094533. 99+100950.1 

SDSS J094533. 99+100950.1 was serendipitously discovered 
by us while searching the full SDSS DR5 (Adelman- 
McCarthy et al. 2007) spectral database with the aim of 
finding interesting quasars for which atypical line and/or 
continuum properties would prohibit them being found and 
classified as quasars by the SDSS automatic pipeline. In par- 
ticular, we looked for quasars with extreme blue continua 
(these would have been missed due to colour cuts used by the 
SDSS to avoid regions of colour space dominated by white 
dwarfs, A-stars and M star+white dwarf pairs - see Richards 
et al. 2002) and quasars with normal (unreddened by dust) 
continua but with non-standard emission line properties. We 
analyzed all, 1 048 960, spectra available in the SDSS DR5 
spectral database, classified by the SDSS automatic proce- 
dures into: galaxies, quasars, stars, sky, and unclassifiable 
objects (see Adelman-McCarthy et al. 2007). 

The spectrum of SDSS J094533.99+100950.1 appeared 
as particularly interesting, showing a normal quasar con- 
tinuum but with atypical, weak emission lines, with practi- 
cally nonexistent C IV and C III] , and weak, although clearly 
visible, Mgll. The object was identified as a quasar by 
the SDSS pipeline due to the presence of the Mgll line, 
and is hence also included in the Schneider et al. (2007) 
SDSS quasar catalogue. The spatial coordinates are a — 
146.39163, 5 = 10.163927 and the determined redshift is 
z = 1.66160 + 0.00025. The quasar's spectrum was observed 
on April 25, 2003 and classified as "excellent" by the SDSS 
database. 

2.1 Emission line measurements - modelling of 
the spectrum 

In order to study the weak emission line properties of our 
WLQ we modeled the spectrum by fitting a power law to the 
underlying continuum, accounting for blended iron emission 
(the iron template of Vestergaard & Wilkes 2001 was used) 
and fitting single Gaussians to the emission and absorption 
lines. Atmospheric emission lines that still remained after 
the automatic pipeline reduction were first removed, and 
the spectrum was corrected for Galactic reddening using the 
Cardelli et al. (1989) extinction curve, adopting a colour 
excess of E{B — V) = 0.06 determined from the Galactic 
neutral hydrogen column from Dickey & Lockman (1990) 
and Stark et al. (1992), and assuming a fixed conversion of 
N{HI)/E{B -V) = 5.0 X lO^^cm"^ mag"^ (Burstein & 
Heiles 1978). 

The underlying power law continuum was fitted to re- 
gions of spectrum uncontaminated by emission lines and 
away from blended iron emission. The continuum windows 
used are shown in Tablf(T] (model "Cont" ) and taken from 
Forster et al. (2001). Blended iron emission was modeled us- 
ing the UV iron template of Vestergaard & Wilkes (2001) 
to spectral windows defined in Forster et al. (2001) (Tabl^ 
model "Fel"). We also experimented with slightly modified 



Table 1. Continuum and iron fitting windows 



Model 


Fitting windows 


Name 


rest frame wavelength range in A 


Conti 


1455-1470, 1690-1700, 2160-2180, 2225-2250, 




3010-3040, 3240-3270 


Fel2 


2020-2120, 2250-2340, 2440-2650, 2900-3000 


Fe2^ 


2020-2120 2250-2300 2500-2650 2850-3000 


FeS* 


1490-1502, 1705-1730, 1760-1800, 2250-2320 




2470-2625, 2675-2715, 2735-2755, 2855-3010 


ContFe4'"'' 


1455-1470, 1685-1700, 2020-2120, 2190-2210 




2250-2300, 2450-2650, 2850-3070 


ContFeS'' 


1440-1470, 1700-1820, 1950-2400, 2450-2700 




2850-3100 


ContFee'^ 


1455-1470, 1690-1700, 2020-2120, 2160-2180 




2225-2650, 2900-3000, 3010-3040, 3240-3270 




3350-3400 



Column (1) - name of model; column (2) - wavelength range of 
spectral windows in which the model was fit. 
^ - continuum fitted as a power law in Forster et al. (2001) 
continuum windows; 

^ - iron template fitted in Forster et al. (2001) iron spectral 
windows, excluding wavelength range 2340-2440A (atmosphere 
artifact -|- feature at 2350A); 

^ - iron template fitted in Forster et al. (2001) iron windows, 
with regions around C II] A2327 and O II A2441 excluded, and 
2850-2900A window (at red side of Mgll) added for a more 
consistent fit; 

* - iron template fitted in spectral windows defined by 
Vestergaard & Wilkes (2001), with regions of CII] emission and 
narrow absorption at CIV (1427-1505A) and Mgll 
(2715-2735A) excluded; 

^ - power law continuum and iron emission fitted simultaneously 
in the "Cont" and "Fel" windows (some windows are narrower 
or slightly shifted to better match the model with the spectrum); 
^ - power law continuum and iron emission fitted simultaneously 
in optimum windows chosen by us and created by comparison of 
the Richards et al. (2003) composites and Vestergaard & Wilkes 
(2001) iron template; 

^ - power law continuum and iron emission fitted simultaneously 
in the "Cont" model windows complemented with a 3350-3400A 
window and iron windows from Forster et al. (2001) . 



iron windows which excluded spectral regions of C II] A2327 
and OIIA2441 emission and/or CIV and Mgll NAL ab- 
sorption (models "Fe2" and "Fe3" in TablfjlJ. We present 
the results of our continuum and iron modelling in Tabled 
In model no. 1 (Cont-Fel-best in Table ^ we fol- 
lowed the procedure described in Vestergaard & Wilkes 
(2001), where first the underlying power law continuum 
was fitted and subtracted from the spectrum, followed by 
the modelling of iron emission. The Vestergaard & Wilkes 
(2001) iron template spectrum was broadened by convolving 
with Gaussian functions of widths between 900 kms^^ and 
9000 kms~^ and separated by steps of 250 kms~^, while 
conserving the total fiux in each template. We chose the 
broadening value for which was lowest. Finally, we sub- 
tracted the best iron fit from the spectrum and again fit- 
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SDSS J094533. 99 + 100950.1 
composite no 3 of Richards et al. 
binned SDSSJ094533 
continuum of SDSS J094533 
iron emission of SDSS J094533 




Figure 1. The restframe spectrum of SDSS J094533. 99+100950.1 corrected for Galactic reddening (yellow line), the binned spectrum of 
the quasar (red line), its fitted, underlying continuum (blue line) and iron emission (green line). For comparison, in black, the Richards 
et al. (2003) composite spectrum (no. 3) is shown, greyshifted to match the SDSS J094533. 99+100950.1 spectrum. 



ted the continuum. This procedure was repeated two more 
times. Final parameters are presented in Table [2] 

In models no. 2-4 (Cont-Fel,2,3) a similar procedure 
was used for the continuum fitting, however broadening of 
the iron template and the scaling factor of the amplitude 
was estimated automatically. We experimented with stan- 
dard and slightly modified iron windows (models Pel, Fe2, 
Fe3 from TabletH). 

In models no. 5-7 (i.e. ContFe4,5,6) power law contin- 
uum and iron template were fitted simultaneously (as a sum) 
to spectral windows including standard and slightly modi- 
fied continuum and iron windows from Forster et al. (2001) 
(models "ContFe4", "ContFe5" and "ContFe6" in Tabled. 

The best continuum and iron emission fit, is represented 
by model no. 1 (Cont-Fel-best). It has the lowest reduced 
value and visually gives the best fit to iron emission in the 
vicinity of the Mg II line, which then leads to the best Mg II 
emission line fit (as shown by the lowest and weakest 
residuals) . In Fig. [1] we plot this continuum and iron emis- 
sion fit over the restframe. Galaxy dereddened spectrum of 
SDSS J094533.99+100950.1. The fit is also included in Fig.[2] 
and the emission line and central engine properties derived 
further in the text are made using this model. 

The best fit continuum has a power law index, a\ 



(fx oc i^"^) equal to —1.232 ± 0.017, which is equiva- 
lent to = -0.77 (/„ oc !/°-). Therefore, the SDSS 
J094533. 99+100950.1 continuum is, in this respect, almost 
identical to the no. 3 quasar composite from Richards et al. 
(2003). However, we also compared quantitatively the under- 
lying continua of SDSS J094533.99+100950.1 and Richards 
et al. (2003) composites no. 1-4 by fitting the latter with a 
power law (to the same spectral windows as for our WLQ 
- the "Cont" model from Table 1) and allowing for inter- 
nal reddening. We checked both Gaskell et al. (2004) and 
SMC-bar (Sofia et al. 2006) extinction curves. We modified 
the composites and allowed for either reddening or dered- 
dening. The lowest was obtained for composite 4 dered- 
dened by E{B - V) = 0.75 for GaskeU et al. (2004) ex- 
tinction, which is relatively high due to its weak depen- 
dence on wavelength, which thus requires large reddening 
to produce any significant change in spectral shape. For the 
case of SMC-bar extinction curve a lower value is needed, 
E{B -V)= 0.02. For composite 3 of Richards et al. (2003) 
no reddening/dereddening was required an acceptable fit. 
For the other two composites reddening was required but 
the fit was always worse than for pure, unreddened compos- 
ite 3. Since the exact wavelength dependence of the intrinsic 
reddening in quasars is still under debate, we decided to use 
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composite 3 for comparison. In Fig. [T] we show the Richards 
et al. (2003) composite no. 3. plotted over the rest frame 
SDSS J094533.99+100950.1 spectrum, corrected for Galac- 
tic extinction. 

Despite the almost identical continua, the two spectra 
are widely different with respect to the emission line proper- 
ties. The Mg II emission line is clearly present but the narrow 
component in the line profile is absent. The contribution of 
the Fell and Felll is almost typical, while the SiIII]H-CIII] 
and C IV emission is extremely weak compared with a typ- 
ical quasar. There are also two deep but narrow absorption 
components (clearly resolved into doublets in the unbinned 
data plot), blueshifted by ~ 8300 km s~^with respect to the 
C IV and Mgll emission line position expected from WLQ's 
adopted redshift. 

The emission and absorption lines were fitted using 
Gaussians and the fits were performed for all continuum and 
iron models defined in Tabled The different models had no 
significant effect on the absorption line fits (resulting in less 
than 11% difference in equivalent widths) and so in Table [3] 
we give the values for one model only Cont-Fel-best. 

The Mgll emission line was more influenced by the 
choice of continuum and iron models, and the equivalent 
width varied from fit to fit by a few percent, with a 15% 
difference between the Cont-Fel-best and the extreme value 
(see Table [4]). The measured kinematic width depended on 
the spectral decomposition, and varied from 5390 km s~^up 
to 6290 km s~^, depending on the model used. However, in 
all studied cases the line was relatively strong and broad. 
Therefore, we consider that the properties of the Mgll line 
are determined reliably which is important from the point 
of view of black hole mass determination. The CIV line is 
consistently faint in all models, although some fits indicate 
the presence of a very weak line while others give values 
consistent with zero intensity. 

For comparison, in Table [5] we give the typical emission 
line parameters measured for QSOs. Single component fits 
and two (narrow and broad) component fits are given when- 
ever possible. We see that the Mgll line intensity is a factor 
of ~ 2 lower than the mean Mgll line intensity of sources 
found in the Large Bright Quasar Survey (LBQS, Forster et 
al. 2001) or observed with the Faint Object Spectrograph 
on the Hubble Space Telescopes (FOS, Kuraszkiewicz et al. 
2002) as well as the value found for the single component 
Mg II fits to the Richards et al. (2003) steep composite. How- 
ever, C IV line (as well as C III]) is over 20 times fainter. The 
weakness of C IV - typically one of the strongest quasar lines 
- qualifies the object as a WLQ. However, the presence of 
Mg II allows for an estimate of the global physical parame- 
ters of our source. 

We also tested the possibility that the redshift determi- 
nation based on the Mg II line is inaccurate and we allowed 
the line centroid and the position of the iron line complexes 
to shift arbitrarily. This did not improve the fit significantly 
although a possibility of such a shift by a few hundred km 
s~^cannot be rejected. 

Since WLQs have been occasionally suggested to be a 
BAL (McDowell et al. 1995), we also calculated the BI in- 
dex (Weymann et al. 1991) in the vicinity of the Mgll line. 
We obtained a value consistent with zero. Also the Al index 
(Hall et al. 2002) for this quasar is zero if we take the un- 
binned spectrum where NAL doublets are clearly resolved 
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Figure 2. The regions of the CIV, SiIII]+CIII] and Mgll lines 
in the spectrum of SDSS J094533.99-I-100950.1 (blue). Green line 
shows the broadened (FWHM(Fe II)=1650 km s~^) iron emission 
template, red and cyan show fitted emission and absorption lines 
respectively, while the black line is sum of all components. Zero 
velocity was at the theoretical wavelength of each line at a given 
redshift. The plotted flux was normalized by the continuum. 
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Table 2. Properties of fitted continuum and iron emission template for different fitting windows. 



No 




continuum 








Fe 11 template 






fitting 






"A 


F 

3000A, 
[erg s cm ' 


A- 




FWHM 

[km s--'-] 


scale factor 


REW 

[A] 




windows 


1. 


Cont-Fel- 

Cont: 
Fe: 


-best 

-1.232 ± 0.017 


(5.24 ± 0.69) ■ 


10- 


16 


1650 ± 130 


0.01176 ±0.00027 


82 


1.51 
1.30 


Cont 
Fel 


2. 


Cont-Fel 

Cont: 

Fe: 


-1.233 ±0.017 


(5.24 ± 0.69) ■ 


10" 


-16 


2150 ± 380 


0.01139 ±0.00027 


79 


1.56 
1.30 


Cont 
Fel 


3. 


Cont-Fe2 

Cont: 

Fe: 


-1.233 ±0.017 


(5.24 ± 0.69) ■ 


10- 


-16 


2680 + 420 


01 208 + 00027 


84 


1.58 
1.30 


Cont 
Fe2 


4. 


Cont-Fe3 

Cont: 

Fe: 


-1.240 ±0.017 


(5.24 ± 0.69) ■ 


10- 


-16 


3440 ± 550 


0.01070 ± 0.00027 


75 


1.58 
1.57 


Cont 
Fe3 


5. 


ContFe4 
Cont-I-Fe: 


-1.194 ±0.011 


(5.40 ± 0.43) ■ 


10- 


-16 


1140 ± 190 


0.00849 ± 0.00043 


58 


1.55 


ContFe4 


6. 


ContFeS 
Cont-I-Fe: 


-1.243 ± 0.010 


(5.22 ± 0.40) ■ 


10- 


16 


2500 ± 300 


0.01195 ±0.00053 


85 


1.71 


ContFe5 


7. 


ContFe6 
Cont-I-Fe: 


-1.215 ± 0.012 


(5.35 ± 0.43) ■ 


10- 


16 


1179 ± 190 


0.01083 ± 0.00043 


74 


1.69 


ContFe6 



Column (1) - model number; column (2) - model name: fits no. 1-4 have continuum modeled as a power law in spectral windows 
defined by model "Cont" in Table 1, independently from the iron emission fitted in iron spectral windows defined by models Fel, 2, 3 in 
Table 1; fits no. 5-7 have continuum and iron emission fit simultaneously to windows defined by models ContFe4,5,6 in Table 1. 
Column (3) lists the power law index, ax (where fx oc A"-^). Column (4) shows the amplitude of the continuum measured at 3000 A. 
Columns (5)-(7) list rest frame full width at half maximum (FWHM) of the Fell template, scaling factor, and Fell rest frame 
equivalent width (REW), respectively. The reduced of the fitting procedure is shown in Column (8). Column (9) lists the fitting 
windows used by us and defined in Table [T] 



Table 3. Rest frame absorption lines properties of SDSS J094533. 99+100950.1 



Line 


[A] 


A(, 

[A] 


REW 

[A] 


FWHM 

[km s" -*-] 


Shift 

[km s-i] 


^abs 


fit: Cont-Fel-best 

MgH 

CIV 


2796.357 
2803.536 
1548.189 
1550.775 


2721.30 ±0.51 
2728.26 ±0.58 
1506.39 ±0.22 
1509.24 ± 0.24 


_o 70 + 0.34 
'='-0.31 

.,c+0.42 
0.00_Q 47 


343 ± 82 
290 ± 110 
455 ± 63 
269 ± 69 


-8268 ± 57 
-8272 ± 64 
-8318 ± 44 
-8251 ±49 


1.59016 
1.59013 
1.58975 
1.59031 



Column (1) - name of absorption line, column (2) - the wavelength of the line seen in laboratory. The absorption lines were fitted with 
Gaussians to a dereddened, rest frame. Foil subtracted quasar spectrum. Column (3) - wavelength of the minimum intensity of the line 
in the fit. Column (4) - Rest frame Equivalent Width (REW) of the doublet (computed as a sum of two Gaussians). Column (5) Full 
Width at Half Maximum (FWHM) of the absorption line. Column (6) lists the line blucshift in km s"^. We present measurements for 
the Cont-Fel-best model from Tabic 2 only, as other models gave similar results. 



- Fig. [2] There is also no sign of a CIV BAL, where the 
observed spectrum allows detection of such a component 
blueshifted up to ~ 20 000 km We conclude that SDSS 
J094533.99+100950.1 is not a BAL or a mini-BAL object. 



archives and software of the Virtual Observatory (VO)0 and 
the Galaxy Evolution Explorer (GALEX) database. This 
included optical u, g, r, i and z photometry from SDSS, and 
near-IR J, H and Kg photometry from the Two Micron All 
Sky Survey (2MASS) and near- (NUV) and far-ultraviolet 



2.2 Spectral energy distribution 



We compiled the broad band spectral energy distribu- 
tion (SED) of SDSS J094533. 99+100950.1 using the data 



^ http://www.euro-vo.org/ 
^ http://galex.stsci.edu/ 
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Table 4. Rest frame emission lines properties of SDSS J094533. 99+100950.1. 



Line 




Ao 

[A] 


REW 

[A] 


FWHM 

[km s-l] 


2cm 


fit: Cont-Fel 

Mgll A2800 
CIII] A1909 + 
C IV A1549 


-best 

Si III] A1892 


2800.29 ± 0.79 
1906.7 ±2.3 
1544.9 ± 2.3 


15.3+}} 

1 40+0.82 


5490 ± 220 
4800 ± 4800 
4000 ± 1100 


1.66187 
1.65845 
1.65451 


fit: Cont-Fel 
Mg II A2800 
CIII] A1909 + 
CIV A1549 


Si III] A1892 


2800.60 ± 0.81 
1905.6 ± 2.4 
1545.6 ± 2.2 


15.5+}'} 
1 69+°-^^ 


5690 ± 230 
5000 ± 5000 
3800 ± 1100 


1.66217 
1.65680 
1.65572 


fit: Cont-Fe2 
Mgll A2800 
CIII] A1909 + 
C IV A1549 


Si III] A1892 


2800.85 ± 0.83 
1904.6 ± 2.5 
1545.9 ± 2.3 


15.21};} 
1 69+°'^* 


5750 ± 240 
5100 ± 5100 
3800 ± 1100 


1.66240 
1.65539 
1.65632 


fit: Cont-FeS 
Mgll A2800 
CIII] A1909 + 
C IV A1549 


Si III] A1892 


2801.64 ± 0.92 
1903.9 ± 2.7 
1545.7 ±2.5 


i5.9i}:^ 


6290 ± 270 
5200 ± 5200 
3600 ± 1200 


1.66316 
1.65455 
1.65597 


fit: ContFe4 
Mgll A2800 
CIII] A1909 + 
CIV A1549 


Si III] A1892 


2799.28 ± 0.91 

1905.6 ± 2.5 

1543.7 ±2.6 


14.0l}-2 
1 07+0.59 
^^•■^'-0.48 
1 50+0.91 


5670 ± 250 
4400 ± 4400 
4300 ± 1300 


1.66091 
1.65680 
1.65252 


fit: ContFc5 
Mgll A2800 
C III] A1909 + 
CIV A1549 


Si III] A1892 


2804.77 ±0.76 
1902.5 ± 2.5 
1547.4 ±2.4 


16.0+}} 
1 01 +0.68 


5600 ± 220 
5400 ± 5400 
3400 ± 1100 


1.66614 
1.65252 
1.65888 


fit: ContFc6 
Mgll A2800 
C III] A1909 + 
C IV A1549 


Si III] A1892 


2798.84 ± 0.87 
1906.0 ±2.4 
1543.7 ±2.6 


14.11};} 

1 oQ+0.58 
J-.oo_Q 47 

1 45+0.89 


5390 ± 240 
4300 ± 4300 
4200 ± 1300 


1.66050 
1.65738 
1.65244 



Column (1) - name of emission line, column (2) - wavelength of the maximum intensity of the line in the fit (emission lines were fitted 
with a Gaussian to a dereddened, rest frame, Fell subtracted quasar spectrum), column (3) - Rest frame Equivalent Width (REW) of 
the line, column (4) - Full Width at Half Maximum (FWHM) of the emission line. Column (5) redshift. 
The fitting window for Mg II is 2750-2840A, for C IV 1530-1580A, and for C HI] + Si III] 1880-1940A. 



(FUV) photometry from GALEX (Edge, private communi- 
cation) . 

No other spectra or photometric data points were found 
in the VO within the matching radius of 5". All available 
images of the quasar and its vicinity were downloaded. The 
quasar is visible in the SDSS and 2MASS images, and also on 
the Palomar Observatory Sky Survey (POSS) plates in the 
optical and near-IR. No X-ray data (XMM-Newton, RXTE, 
INTEGRAL) were found. We searched for X-ray spectra in 
the ROSAT catalogues, but found no sources with a flux 
significantly higher than the X-ray background (Soltan, pri- 
vate communication). The radio properties of our quasar and 
its neighbourhood were checked in the Parkes-MIT-NRAO 
(4850 MHz survey), and in the Faint Images of the Radio 
Sky at Twenty-centimeters (VLA FIRST, 1.4 GHz). In both 
cases the radio background was found and the spatial reso- 
lution was too low to observe the quasar. The same results 
were obtained for IRAS and EUVE catalogues - the signal 
did not differ significantly from the background. 



The SED of SDSS J094533.99+100950.1 is shown in 
Fig. |3] It does not differ significantly from typical quasar 
SEDs (e.g. Elvis et al. 1994, Richards et al. 2006), consistent 
with the conclusion reached by Diamond-Stanic et al. (2009) 
for other WLQs. The deviation from the Elvis et al. (1994) 
median SED in the far UV (at NUV and FUV GALEX data 
points) is due to the simple interpolation between the UV 
and soft X-rays in the Elvis et al. (1994) SED, which does 
not constrain the true QSO SED at A < 1000 A. However, 
lower flux in the far-UV remains a possibility in our WLQ, 
which may lead to low CIV emission due to lower ionizing 
flux. 



3 GLOBAL PARAMETERS OF SDSS 
J094533. 99+100950.1 

Our WLQ has a monochromatic luminosity at 3000A of 
ALa(3000A) = 2.95 x 10^*^ erg s"\ (assuming Ho = 71 km 
s"^ Mpc"\ SIa = 0.73 and fi^ = 0.27 - Spergel et al. 2007), 
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Table 5. Emission line properties in other quasar samples. 



Line 


REW(composite) 

[A] 


REW(FOS) 

[A] 


FWHM(FOS) 

km s^^ 


REW(LBQS) 

[A] 


FWHM(LBQS) 

km s~^ 


Mgll 
single 
broad 
narrow 


30.14 


136 ±92^ 
33 ±8 
29 ±9 


3840 ± 1850 
8230 ± 1850 
3260 ± 760 


39 ±22 
37± 22 
28± 13 


5160 ±120 
8660 ± 300 
3510 ±110 


Cm] 


21.45 


21 ± 5 


4890 ± 780 


28 ± 15 


7820 ± 170 


Civ 
single 
broad 
narrow 


26.86 


21 ± 16^ 
63 ±7 
30 ±4 


4430 ± 3070 
11090 ± 1070 
2900 ± 290 


38 ±20 
45 ±23 
17.8 ±9.3 


7720 ± 150 
10960 ± 360 
2860 ±110 



Emission lino parameters of SDSS J094533. 9+100950 are compared with: in column (2) the Richards et al. (2003) composite no. 3 
(Fig. [1} built on the basis of 770 quasar spectra, in columns (3), (4) the 158 AGN observed with Faint Object Spectrograph on the 
Hubble Space Telescopes (FOS; heterogeneous sample of mixed type AGN, Kuraszkiewicz et al. 2002) and in columns (5), (6) the 993 
quasars from the Large Bright Quasar Survey (LBQS; Forster et al. 2001). In Column (1) "single" means that the line was fitted by a 
single Gaussian, while "broad" and "narrow" refer to the broad and narrow components of the line. 

^ - The single Gaussian fits to Mg II and C IV emission linos in FOS were based only on a few, low quality spectra and thus may not be 
representative. 
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Figure 3. The broad band SED of SDSS J094533.99+100950.1. 
The photometric IR data point (triangles) come from 2MASS 
(J,H,Ks colours). Optical photometry (squares) is from SDSS 
(u,g,r,i,z filters). The NUV and FUV points (hexagons) are from 
GALEX. The spatial coordinates of the source differ between 
SDSS and 2MASS by ~ 0.39 arcsec and between SDSS and 
GALEX by ~ 0.48 arcsec. We assume the data from SDSS, 
2MASS and GALEX come from the same object. Galactic dered- 
dening was done using the values E(B-V)=0.060, Ry = 3.2 and 
formulae of Cardelli et al. (1989). Neither intrinsic dereddening of 
the quasar nor subtraction of Fe II were done. For comparison the 
line shows the SED from Elvis et al. (1994) marking with dashed 
line the frequency range of interpolation between the far UV and 
soft X-ray band. 



which is high but not extremely high for a quasar. Using a 
bolometric correction BC3000 ~ 5.15 determined at 3000 A 
(Shen et al. 2008 and Labita et al. 2009, based on Richards 
et al. 2006) we obtain a bolometric luminosity of 1.52 x 10*^ 
erg s~^. 

We use the following equations from Kong et al. (2006) : 
to determine the central black hole mass: 



Mbh = 3.4x10 



Mbh = 2.9x10 



XL 



(lO^^ergs-i 



FWHMMgii 



Y L 



Mgll 



\ lO'^^ergs" 



1000 kms- 



FWHMmsII 



1000 kms-i 



M0(1) 



Ms (2) 



where XL <■ and LmeII are the monochromatic continuum 

3000A ^ 

luminosity at 3000 A and the Mgll line intensity, respec- 
tively. Using the Mgll width from the best fit (Cont-Fel- 
best in Tabled]) we obtain a black hole mass of 2.7 x l(f Mq 
from equation This is a typical value for a distant quasar 
(e.g. Kelly et al. 2008, Vestergaard et al. 2008). Using equa- 
tion ([2| we obtain a slightly lower value of the black hole 
mass, of 1.5 x 10® Mq. Similar values of Mbh are obtained 
when using formulae from McLure & Dunlop (2004). 

The first value of the black hole mass implies the Ed- 
dington ratio of 0.45. This value is very conservative from the 
point of view of possible super-Eddington accretion, which 
is possibly present in some AGN (Collin & Kawaguchi 2004, 
Jin et al. 2009). It is consistent with several other determi- 
nations of this ratio in quasars. Quasars in the redshift range 
between 1 and 2 in the study of Sulentic et al. (2006) have 
Eddington ratios between ~ 0.2 and ^ 1, the value 0.25 is 
the mean favoured by statistical studies of Shankar et al. 
(2008). Labita et al. (2009) found a value of 0.45 as a statis- 
tical upper limit of accretion rate in their sample of quasars. 
The second value based on black hole mass from equation ((2| 
gives a somewhat higher ratio of 0.79, which is less reliable 
because of the weakness of the Mg II emission line. We also 
computed the mean and standard deviation of the Edding- 
ton ratio for all models in Tableland obtained the following 
means: 0.429 ± 0.037 using equation (O and 0.752 ± 0.077 
for equation ((2|. In the newest studies of the relation be- 
tween Mgll line parameters and black hole mass - Wang 
et al. (2009) do not use Mgll equivalent width but rather 
the Mgll FWHM and continuum luminosity at 3000A. Us- 
ing the relation from Wang et al. (2009) the mean black hole 
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mass of SDSS J094533.99+100950.1 is (3.09±0.23) x 10^ Af© 
and the mean Eddington ratio is 0.396 ± 0.026. 

Additionally, since the Mgll emission line intensity is 
lower than in typical quasars, the second method is probably 
less reliable than the first, based on the continuum, and so 
the Eddington ratio of 0.45 is more likely. 



4 DISCUSSION 

The weak line quasar SDSS J094533. 99+100950.1 is an ex- 
ceptional object with almost nonexistent C IV emission, nor- 
mal Fell emission and visible, although somewhat fainter 
and atypical, Mg II emission line. The presence of iron emis- 
sion with standard intensity argues against any relativistic 
enhancement of the continuum in this source. The quasar 
continuum is normal, similar to the steep composite no. 3 of 
Richards et al. (2003). 

The lack of strong emission lines (except for Mgll) is 
likely to be intrinsic, consistent with the conclusion reached 
by Diamond-Stanic et al. (2009) for other WLQs. A remark- 
able property of the SDSS J094533.99-f 100950.1 spectrum 
is the presence of Low Ionization Lines (LILs) such as Mgll, 
which are thought to form close to the accretion disk surface 
(see e.g. CoUin-Souffrin et al. 1988) and the absence of High 
Ionization Lines (HILs) such as C IV. The narrow compo- 
nent of Mgll, which is thought to form at larger distances 
from the nucleus, is also absent from the spectrum. 

McDowell et al. (1995) listed 10 possible explana- 
tions for the weakness of emission lines in their WLQ, 
PG 1407-1-265, but rejected most of them on the basis 
of overall quasar properties. Among those remaining are: 
high extinction or BAL effects, which are inconsistent with 
SDSS J094533. 99+100950.1 properties, and abnormal ioniz- 
ing continuum, which remains a possibility. In particular, the 
option of high (super-Eddington) accretion rate may seem 
attractive since the determination of the Eddington ratio 
for our object allows (but does not require) for such a pos- 
sibility (see Sect. (3)1. However, the expected trends found 
by PGA analysis performed on the optical/UV emission line 
properties of PG QSO (Shang et al. 2003; see Fig. 10) show 
that, with increasing Eddington ratio, the UV Fell emission 
(measured around Mg H) decrease^. Si HI] +C HI] emission 
increases, and the widths of Mgll and SiIII]+GIII] lines 
decrease. The opposite is seen in our object. Additionally 
the Si IV line should become stronger, but this is not seen 
in the WLQs of Shemmer et al. (2009). Since the emission 
line properties of WLQs do not follow the trends for high 
L/LEdd objects predicted by PGA analysis, something dif- 
ferent than high L/LEdd is causing the atypical emission line 
behaviour. 

The only explanation, which in our opinion is currently 
consistent with the intriguing emission line properties of our 
WLQ is that the quasar activity in this source has just be- 
gun. 

The formation of an accretion disk in an AGN is a long 
time phenomenon, occurring in millions of years (the vis- 
cous timescale in the outer disk, Siemiginowska et al. 1996) . 

^ not increases, as is the case for optical Fell measured around 



However, when the disk finally approaches a black hole, i.e. 
when the inner radius moves from ~ 10 — 20-Rschw to ISGO 
(Inner Stable Gircular Orbit), the disk's luminosity rapidly 
increases on a hundred year timescale and X-ray emission 
starts (Gzerny 2006). 

The immediate stage after is the irradiation of the outer 
disk which happens in the light travel time across the re- 
gion (years). At this stage the accretion disk continuum ap- 
pears already to be typical for an AGN. However, neither the 
broad-line region (BLR) nor the narrow-line region (NLR) 
region, exist. 

If our explanation is correct, our WLQ should have no 
narrow [O HI] A5007 emission but its LIL H/3 line should have 
already formed (as has the Mgll line). For our z = 1.6 ob- 
ject these lines lie outside the observed wavelength range 
of the optical telescopes used by SDSS, but could, in prin- 
ciple, be observed using an IR facility. In another WLQ, 
PG 1407+265, (McDowell et al. 1995) the H^ line is weak, 
contrary to our expectations. However, significant influence 
of the atmosphere above 9000A leading to low spectrum 
quality at those wavelengths may partially hide the line (see 
Fig. 1 of McDowell et al. 1995). The puzzling properties of 
the X-ray emission of PG 1407+265 may be consistent with 
the young age of the source as the emission significantly 
varies in time, and is only occasionally dominated by a jet 
(Gallo 2006) indicating that a stable jet has not yet been 
launched. 

A relatively young age was also suggested for the class 
of Narrow Line Seyfert 1 (NLSl) galaxies, in the context of 
cosmological evolution (Mathur 2000) but their BLR is well 
developed. WLQs thus possibly represent an even earlier 
stage but not necessarily in the sense of a global evolution 
since the argument of their young age is based on BLR for- 
mation stage. If the observed activity is the first active stage 
of the nucleus then indeed the WLQ stage should be seen 
first, later replaced by the NLSl stage (with efhcient growth 
of the central black hole), followed by the broad emission line 
stage (quasar or Seyfert 1). However, it is also possible that 
the observed active stage represents the source's reactiva- 
tion. The amount of intrinsic (nuclear) dust along the line 
of sight does not seem to be large and the past activity may 
be the natural explanation for such a clean environment for 
an otherwise young object. The large mass of the black hole, 
particularly in PG 1407+265 (logAfen = 9.8 based on CIV 
FWHM and the formula of Vestergaard & Peterson 2006) 
might rather indicate a relatively late evolutionary stage, 
favouring the reactivation scenario. A few stages of quasar 
life are also predicted by Hopkins & Hernquist (2009). They 
postulate that the life of quasars consists of a merger and a 
few episodes of high accretion rate. The activity is likely to 
cease through a "naked" AGN stage (Williams et al. 1999, 
Hawkins 2004) since when the accretion rate drops, the op- 
tical disk recedes outward and the component of the BLR 
linked to the disk wind disappears (e.g. Gzerny et al. 2004, 
Elitzur & Ho 2009). 

An intermittent character for the quasar activity is also 
more consistent with the statistics of WLQ. The sources 
are rare but not as rare as we might expect (based on the 
following discussion). The irradiation causes the formation 
of a wind in the upper disk atmosphere, and this material 
slowly rises up in the vertical direction. This wind is cus- 
tomarily now considered as the source of the material from 
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BLR clouds (Murray et al. 1995; Proga et al. 2000; Che- 
louche & Netzer 2003; Everett 2005). The initial rise is not 
highly supersonic. Assuming velocities ~ 100 km s~^ (e.g. 
Proga 2003, Everett 2005) we can estimate how much time 
the material takes to depart significantly from the disk's 
surface layers. Assuming a disk radius of 1 pc, a rise by ten 
percent of the radial distance takes 1000 years, and at that 
stage the low ionization lines forming close to the disk, like 
Mgll and H/? appear. The rise to the height comparable to 
the radial distance takes a factor 10 times longer, and only 
after that time the highly ionized lines such as C IV ap- 
pear in the spectrum. The narrow emission line components 
take longer to form and are likely to be absent at the ini- 
tial stages. A similar argument was used for the explanation 
of the relative faintness of the [O III] line in GPS quasars 
which are also thought to be relatively young (Vink et al. 
2006). Thus, a WLQ phase lasts for about 1000 years. There 
are almost 100 quasars classified as WLQ (Shemmer et al. 
2009, Diamond-Stanic et al. 2009, Plotkin et al. 2010) in the 
whole sample of ^ 100 000 SDSS quasars. However only a 
fraction of those, including SDSS J094533. 99-1-100950.1 (this 
paper), PG 1407+265 (McDowell et al. 1995) and HE 0141- 
3932 (Reimers et al. 2005), may represent a more advanced 
stage with partially developed low ionization lines. There- 
fore, the probability of observing this evolutionary stage can 
be estimated as 10""^ — 10~*. This would indicate that the 
typical duration of a quasar's active phase is only ~ 10^ — lO'' 
years. On the other hand, the lifetime of the luminous phase 
of quasars is ~ 10^ < fqso <~ 10* yrs (Haiman & Hui, 
2001), but in some sources extends over a period ~ lO'' yrs 
(Martini & Weinberg, 2001). The two numbers can be recon- 
ciled if the total active phase of a quasar consists typically of 
100 separate subphases, each of those starting with a WLQ 
stage of BLR buildup. 

The increasing probability of occurrence of WLQ among 
high redshift objects (Diamond-Stanic et al. 2009) agrees 
with our hypothesis if, in the earlier Universe, minor mergers 
are more frequent and the typical single active subphase 
lasts shorter. 

Thus the SDSS J094533.99+100950.1 quasar with its 
remarkable properties may have a Rosetta stone impact on 
our understanding of quasar activity cycle and the nature 
of WLQ. 

Finally, it is important to note that the definition of 
a WLQ has not been yet worked out. Some objects with 
extremely faint/nonexistent lines, may indeed be BL Lacs. 
Others, like PHL 1811 (Leighly et al. 2007), considered as 
a WLQ prototype, may be explained by higher Eddington 
ratio since the lines are not only fainter but also narrower 
(Mgll line of 2550 ± 110 km s"^ vs. 5490 ± 220 km s"^ 
in J094533.99-f 100950.1). Our hypothesis of young/newly 
active sources would not apply to those classes of objects. 



5 CONCLUSIONS 

Quasar SDSS J094533.99-K100950.1 is an exceptional exam- 
ple of the Weak Line Quasar class. The broad band spec- 
trum is typical for quasars, very similar to the second steep- 
est composite spectrum of Richards et al. (2003). It has a 
relatively well developed broad Mgll line and the typical 
Fell contribution but its high ionization lines, like CIV, are 



weak/absent. The BAL phenomenon is not visible in the 
spectrum so the lack of C IV line is an intrinsic property of 
the source. The presence of Mg II and typical Fe II emission 
excludes a non-thermal, blazar-like origin of the spectrum. 
The significant width of Mgll line (FWHM of 5500 km s"^) 
and the Fell properties argue in turn against an extremely 
high accretion rate in this object. Objects like these are diffi- 
cult to find since a typical search for WLQs, based on finding 
spectra with extremely weak or non-existent emission lines 
is likely to miss sources where one of the lines is relatively 
strong. However, such objects are important to find/study 
as they are likely to shed light on the the mechanism driving 
all (or at least some) WLQs. 

In the present paper we propose the following hypoth- 
esis for the nature of WLQ. AGN are generally classified 
according to their viewing angle and Eddington ratio. How- 
ever, those objects evolve and some evolutionary stages may 
be too short to produce a consistent equilibrium between 
various emitting components in the nucleus like the BLR, 
NLR and the accretion disk. It is an interesting possibil- 
ity that the WLQ class represents such a stage, at the 
onset of quasar activity. The relatively high occurrence of 
these sources may additionally suggest that a quasar's ac- 
tive phase has an intermittent character. The whole active 
phase consists of several subphases, each starting with a slow 
development of the BLR region which manifests itself as a 
WLQ phenomenon. Further observational tests of this hy- 
pothesis are clearly needed. 
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